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Abstract 

An exact inversion formula for the Lorentz integral transform (LIT) is provided together with the 
spectrum of the LIT kernel. The exponential increase of the inverse Fourier transform of the LIT 
kernel entering the inversion formula explains the ill-posed character of the LIT approach. Also 
the continuous spectrum of the LIT kernel, which approaches zero points necessarily to the same 
defect. A possible cure is discussed and numerically illustrated. 
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I. INTRODUCTION 



Integral transformations 

/oo 
dxn(y, x)R(x) (1) 
-oo 

are widespread used in physics. According to Hadamard jj, the problem of finding the 
unknown function R for a given transform L is called "well-posed" , if the following conditions 
are fulfilled: 

• The function R exists and is unique. 

• The function R depends continously on the input L. 

Otherwise, the problem (CQ) is called "ill-posed" . As example, let us consider the well-known 
Laplace transformation 

/•OO 

L(y) = / dxe- xy R{x) (2) 



o 

which is widely used in physics and engineering. It turns out to be ill-posed. For illustration, 
consider the test function 

*M = ^ (3) 

e 

and its Laplace transform 

L(y) = earctan( — ) . (4) 

For e — > 0, the maximum norm of L goes to 0, whereas the maximum norm of R is 1 for 
any e. 

In this paper we study integral transforms, where the kernel k in ([!]) has the formQ] 



n(y,x) = K(z = y - x) (5) 
In detail, we will concentrate on the Lorentz integral transform (LIT) (<7j > fixed) 

Klit(z) = --^- 2 (6) 

which has recently been applied extensively in photo- and electronuclear physics Q]. 

In the LIT-approach, due to the choice Cj ^ 0, the calculation of physical observables 
in the A-particle scattering problem can be traced back rigorously to the solution of an 



appropriate bound state problem. This is of course a tremendous technical simplification 
and offers a unique possibility to carry out rigorous ab initio calculations on few-nucleon 
systems even beyond mass number A = 4 [2]. However, the LIT is known to be ill-posed as 
will be again outlined below. 

In general, in order to obtain R from a given ill-posed transform, an appropriate regu- 
larization scheme is essential. For example, in one class of inversion schemes for the LIT, 
the numerically gained Lorentz transform is expressed as a finite sum of appropriate basis 
functions whose inversion in explicitly known, i.e. 

N 

L(y) = 5> x ,(y) ■ (7) 

1=1 

Due to the ill-posed character of the LIT, the upper value N of basis functions must be 
limited for obtaining reliable results. Otherwise, the gained solution for R may contain 
strong, unphysical oscillations. This reduction in the numerical resolution is nothing else as 
a regularization, see |3|] for further details. 

In this paper, we intend to tackle the inversion problem of any ill-posed integral transform 
of type (j3j) in a conceptually completely new manner, namely, without any use of regular- 
ization techniques. For that purpose, we present in section 2 a new inversion formula for 
kernels of the type (jHJ) which will then be applied to the LIT. It turns out that this inversion 
formula directly exhibits the ill-posed character of the LIT. We further provide the spectrum 
of that LIT kernel, which also makes its ill-posed character evident. 

Section 3 is devoted to a numerical case study for the LIT approach, which illustrates 
the very source for the ill-posed property. A possible reduction of that illness is proposed, 
i.e. a strategy will be developed to heal, at least under certain circumstances, the ill-posed 
character of the LIT. We end with a brief summary in section 4. 

II. MATHEMATICAL PROPERTIES OF KERNELS INCLUDING THE LIT 

For kernels of the type fl3j), eq. (pQ) has the form 

L(y) = / dxKiy - x)R(x) (8) 
Jo 

which Fourier transformed leads to 

L(k) = \/2nK(k)R(k) (9) 



with 



/(*) = -= / dye iky f(y) (10) 

V Z7T 



OO 



for any function /. Immediately, from (Q, it follows the inversion in closed form 

Provided that the Fourier transform K(k) of the kernel is a continuous function and that 
the infinum of its absolute value has a lower bound larger zero, i.e. 

C = inf\K\>0, (12) 

we obtain for the L 2 norm 

poo 1 

\\R\\ 2 = dx\R(x)\ 2 < — — ||L|| 2 . (13) 



4tt 2 C 2 

In consequence, we can formulate a very simple criterion to distinguish between ill-posed and 
well-posed integral kernels of the type (jSJ). If C > is fulfilled, the integral transformation 
is well-posed. Otherwise, it turns out to be ill-posed, as it becomes obvious by studying the 
spectrum of integral kernels of the type (j5J). It is defined as 



dxK(y - x)xk (x) = v(k )Xk (y) , (14) 

— oo 

or Fourier transformed as 

V2^K(k)x k0 (k) = v(k )Xk (k) . (15) 

This leads to 

Xk (k)~5(k-k ) (16) 

and 

Ufa) = V2^K(ko) . (17) 

Therefore, the spectrum of any arbitrary integral kernel of the type (jSJ) is given by its 
Fourier transform. If a kernel fulfills C = in (1121) . the spectrum is continuous and has zero 
as accumulation point. In consequence, the kernel is obviously ill-posed and an unprotected 
inversion cannot work. 
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In case of the LIT, K(k) is given by 

K LIT (k) = ^e~^ . (18) 

It fulfills C = 0, its spectrum is continuous, 

ULiAko) = , (19) 

and approaches zero for large ko- Therefore, the LIT is ill-posed. Its eigen functions are 
given by 

Xko (x)~e^ x . (20) 

III. A NUMERICAL CASE STUDY FOR THE LIT 

The application of the inversion formula (ITTj) requires a reliable numerical treatment of 
Fourier transforms. In this work, this is performed with the help of Filon's integration 
formula 4j. In order to check our numerical routines, we take a simple analytical test case 
for the function R: 

R(x) = y/xe- ax Q(x) , (21) 

with denoting the Heavyside step function, i.e. the threshold of this function is placed at 
x — 0. a is a free parameter which we choose from now on as a = 0.05 MeV" 1 . Following the 
nomenclature in literature, we call R from now the response function. Its Fourier transform 
is known analytically 

^»-V / l(^TF) ie " |ata,l( " ) (22) 

so that the quality of our numerical routines for the Fourier transform can be tested 
straightforwardly. For that purpose, we cut at first the integrations bounds in (flUl) according 
to 



R{k) = -= \ dxe lkx R[x) — ► -= / dxe tkx R{x) (23) 

v27T J— oo v27T J—Xr, 



<- m a x 



with x max = 3000 MeV. With respect to the strong exponential decrease of our test 
function (I2ip . this upper value for x max is more than sufficient. The remaining integral 
in (!23j) is now evaluated with the help of Filon's integration formula using iV equidistant 
grid points. It turns out that the numerical precision in calculating the Fourier transform 



increases naturally with increasing N, but decreases strongly with increasing argument k of 
the Fourier transform. For example, for iV = 20001 mesh points the numerical error in the 
real part of the Fourier transform is of the order of about 1.7 percent for k = 0.8 Me V -1 , 
increasing to about 6.9 percent for k = 2.0 MeV -1 . 

At next, let us check the inversion formula (TIT]) for our test response (12 ip . For that 
purpose, we calculate at first numerically the LIT transform L(y) of R(x) and then its 
Fourier transformations L(k) corresponding to (123]) with x max = 3000 and iV = 20001 
equidistant mesh points. Then, with the help of the known Fourier transform of the Lorentz 
kernel ( fl8|) . we can perform the inversion ( TTTT) where we again, similar as in ( 1231) . have to 
cut the integration at the lower and upper limit, i.e. 

j-y, \ If 00 L(k) _ ikx 1 l(k) _ ikx , n ,. 

R ( x ) = 7T / dk ~ ' e — ► — / dk _ . . e . (24) 
2n J-oo K LIT (k) 27r J~k max K LIT (k) 

This second Fourier integral is calculated with Filon, considering again iV = 20001 equidis- 
tant mesh points in the interval [— k max , k max \. 

In order to obtain in general a quantitative estimate for arising numerical uncertainties, 
we define in this context the quantity 

= \\R num (x) - R(x)\\ 2 

\\R num (x)\\ 2 { ' 

with R num denoting our numerical result for the real part of the response derived via the 
inversion formula (fTTj) . Concerning the norm ||, we use f fl3|) . i.e. the L 2 norm. 

In our first attempt, calculation A, we choose k max = k^ ax = 0.8 MeV -1 in (|24l) . The 
resulting response (e = 2.2 10~ 3 , dashed curve in Fig. [TJ, using ff24l) . turns out to be qual- 
itatively correct beyond the threshold region, but unsatisfactory at threshold. The reason 
for this failure and a convenient improvement will be discussed below. Apart from that 
point, we can conclude that our new inversion formula (ITTj) works, at least in principle for a 
LIT with arbitrarily high precision. In practical applications, however, the LIT L(y) is not 
known perfectly well. In order to take this important fact into account in the present case 
study, we proceed as follows. We take, for a given value of y, the exact Lorentz transform of 
( |2~TT) ("exact" within the limits of the compiler precision) in the form . X\X2 . . io- d 

with Xi, d integer numbers and x\ 7^ 0. Then, we substitute this number by hand according 
to 

Q.XiX 2 X' i x i x 5 x & ... 10 _d — ► O.X1X2X3 10~ d , (26) 
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i.e. we cut all relevant digits beyond the first three ones. This ad hoc procedure simulates 
a numerical error in L in the range ~ 1CT 3 — 10~ 2 - a level which can, under favourable 
conditions, reached in state-of-the-art calculations in nuclear physics. This modified Lorentz 
transformation, used from now on in this work, can hardly be distinguished from the original 
one. However, the resulting response function, yielded by the application of fill I) , shows large 
oscillations (e = 5 10~ 2 , dotted curve in Fig. [p. This calculation is denoted from now on as 
calculation B. This numerical fact shows insistently the well-known ill-posed character of 
the Lorentz integral transformation: a small change in L (lhs. of (|T])), generated in our test 
case by the modification ( 1261) . leads to large changes in the response R (rhs. of equation 
©)• 

As a novel feature, our inversion procedure allows now to pin down precisely the source 
for the ill-posed character of the Lorentz transform. For that purpose, Fig. [2] shows for 
calculation A and B the real part of the corresponding Fourier transforms R(k) yielded via 
®, i.e. 



One obtains a remarkable fact. For arguments k > 0.4 MeV -1 , the Fourier transform yielded 
in calculation B (dashed) differs considerably from the "exact" result A (full). This can be 
easily understood by the explicit form of Kurik) ( fT8l) . Its inverse is growing exponentially 
so that in fl27|) even small errors in L(k) are tremendously amplified. This fact also explains 
intuitively why for small arguments k both Fourier transforms are almost identical. This 
important result can also be understood from a more rigorous point of view. For k — > 0, we 
have 

lim V2^R{k) = f dxR(x) = [ dyL{y) , (28) 

where the last equation follows from the structure of Lorentz kernel ([S]). In consequence, if 
the LIT L contains only a small error, then the last integral in fl28|) will also contain only a 
small error. Therefore, the Fourier transform R(k) can be determined with great precision 
for k — > via (Q - despite that the LIT is ill-posed. 

Let us summarize this important result again: The LIT is an ill-posed integral transform. 
However, this ill-posed character does not affect at all the numerical stability of the obtained 
Fourier transform R(k) of the resulting response for sufficiently small arguments k. Only 
for moderate and large arguments k, numerical instabilities arise which can be traced back 
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to the ill-posed character of the LIT. 

In consequence: If it were possible to determine, from general principles, the asymptotic 
behaviour of the Fourier transform R(k), one could hope to circumvent, or at least to reduce 
significantly, any arising numerical instability problems due to the ill-posed character of the 
LIT. From a certain point of view, the ill-posed character of the LIT could be "healed" . 

For our case study, we proceed now as follows. At first, we have to fix the value kthres up 
to which the ill-posed character is irrelevant for the determination of R(k) via ff27|) . Unlike 
as in our case study, the exact Fourier transform R(k) is of course not known in general. 
In practice, one could help oneself by repeating the calculation for a slightly different value 
0"j. In an ideal calculation, without any numerical errors, the resulting Fourier transform 
R(k) should of course be independent from the chosen value crj. In Fig. [2] the real part 
of R(k) can be compared for our standard value Oi = 10 MeV (dashed) to the one for 
(jj = 11 MeV (dotted). Till about k = 0.25 MeV -1 , both Fourier transforms are almost 
identical, however, especially for arguments larger than about k = 0.4 MeV -1 substantial 
differences arise. Therefore, we could trust our results only till about kthres = 0.25 MeV -1 . 
For k > 0.25 MeV -1 we have to substitute the unphysical result obtained in calculation B 
by a more reasonable one. For thatpurpose, we exploit general theorems on the asymptotic 
behaviour of Fourier transforms 3, E|- It turns out that quite in general the latter is 
unambigiously fixed by the knowledge of the threshold behaviour of the response function 
R. Now, the threshold behaviour of R is, for a given reaction, in general known from basic 
physical principles. Without any considerable loss of generality, we can therefore assume 
that we know not only in our case study, but quite in general the threshold behaviour of our 
unknown response (121]) . i.e. 

R(x) ~ Vx . (29) 



Due to 



6|, its Fourier transform must therefore behave asymptotically as 

R(k) ~ -^3 , (30) 



in agreement with the exact analytical result 

Now, it is almost straightforward how to proceed. Only for k < k t hres we use in (El 
numerical results for L(k). For k > k t hres, we use instead the asymptotic ansatz (compare 
with ©) 

L asy {k) = V2^K(k)R asy (k) (31) 
8 



with 

Rasy(k) = . (32) 

The complex parameter c is choosen to guarantee a continous integrand at k — ±k t h res in 
f T2"4"]) . By this prescription, it is of course easily possible to extend the limit k max in ([23]) to 
arbitrarily large values without any arising numerical errors. In the resulting calculation C, 
we choose k max = 10 MeV -1 with again N = 20001 equidistant mesh points in the region 
[~k max , kmax]. The resulting response is depicted as the dashed curve in Fig.[3](e = 4.2 10~ 4 ). 
One readily realizes two very important improvements: (i) The heavy oscillations found in 
calculation B have vanished, (ii) the threshold behaviour is considerably improved compared 
to calculation A. Both features can be traced back to the improved treatment of L(k) by 
exploiting the known asymptotic behaviour of R(k). Since we assume that we know the 
threshold behaviour of the response, the still arising problem of a small, but nonvanishing 
numerical result for R(x) below threshold is irrelevant in practice. 

In general, for an unknown response, it is of course not clear whether the asymptotic 
form is already working for k = k t hres where the Fourier transforms R(k) start to diverge 
from each other for different choices of o"j. The present case, considers, from that point of 
view, even an "almost" worst case scenario, since the Fourier transform is decreasing rather 
slowly, R{k) ~ It turns out that for k = 0.8 MeV -1 , the asymptotic form and the exact 
form of the real part of the Fourier transform, (132]) vs. (122]) . differ by not less than about 
40 percent. Nevertheless, even this quite reduced precision in the knowledge of R{k) for k 
at quite moderate values beyond k t hres is absolutely sufficient to obtain reasonable results 
in the inversion of the LIT. If the Fourier transform decreases stronger, e.g. exponentially, 
one can expect even more reliable results within our proposed inversion technique. 



IV. SUMMARY 



We illuminated the ill-posed character of the LIT approach by providing an exact inver- 
sion formula and giving the spectrum of the LIT kernel. The exact inversion formula requires 
the Fourier transform of the LIT transformed response, which has to be obtained numeri- 
cally, whereas the Fourier transform of the LIT kernel is analytically known. Since K^ 1 {k) 
occurs in the inversion formula, even small numerical errors in L(k) at large /s-values are 
drastically enhanced which lead to violent oscillations of R(x) in the inverse Fourier trans- 
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formation. We illustrated that situation in a numerical test case. Evaluating L(y) and its 
FT very precisely and using the inversion formula we recovered the model response R(x) 
with satisfactory precision, at least beyond the threshold region. However, truncating L(y) 
to about three digits revealed very drastically the ill-posed character of LIT, namely leading 
to wild oscillations in R(x) evaluated by the inversion formula. It is not the behaviour of 
L(k) for small A;- values which causes the oscillations but the behaviour for the large k- values. 
This is clearly demonstrated by putting L(k) to its known asymptotic behaviour for k larger 
a certain k t hres- The resulting R(x) gained through the inversion formula is close to the 
exact one, without oscillations. In this framework, also the threshold behaviour of R(x) can 
be reproduced satisfactorily. 

As not oulined in this work, we verified in addition that for a more structured model 



response with two peaks, like the one used in 



3], that inversion formula works equally well, 



provided one corrects again the asymptotic behavior of L(k). 

In the applications of LIT to electro weak processes in nuclear physics, L(y) can be 
determined in various manners, see [2] for further details. Using our inversion formula, the 
Fourier transformed L(k) of a numerically gained LIT is required. We conjecture that our 
results lead to the requirement to determine L(y) with such a high precision that its Fourier 
transform is well under control for small and intermediate arguments, but not necessarily 
for large arguments. If this is guaranteed, and if the threshold behaviour of the response is 
known, the proposed method allows to heal the ill-posed character of the LIT without any 
use of regularization techniques. From a conceptual point of view, this a very remarkably 
result. 

Since our proposed ansatz is quite general, provided the integral kernel is given in the 
form (jSJ), one can apply the techniques developed in this work to a variety of different other 
ill-posed problems in physics and engineering. 



[1] J. Hadamard, Princeton University Bulletin 13 (1902), 49 

[2] V.D. Efros, W. Leidemann, G. Orlandini, N. Barnea, J. Phys. G: Nucl. Part. Phys. 34 (2007), 
R459 

[3] D. Andreasi, W. Leidemann, Ch. Reiss and M. Schwamb, Eur. Phys. J. A 24 (2005), 361 



10 



2.5 




-0.5 1 1 1 1 1 1 1 

20 40 60 80 

x [MeV] 

FIG. 1: Full curve: Exact Response function (|2Tj) . dashed curve: Response function yielded by 
inversion formula (jlip (calculation A) with k max = 0.8 MeV- 1 . Dotted curve: result for response 
in calculation B with k max = 0.8 MeV -1 . In all curves, <7j = 10 MeV has been used. 
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FIG. 3: Full curve: exact response function (f2Tj) . dashed curve: response yielded in calculation C 
with <Ti=10 MeV. 
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